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We explore patterns of effective restoration of the chiral Uy 4 (l) symmetry using an extended 
three-flavor Nambu-Jona-Lasinio model that incorporates explicitly the axial anomaly through the 
’t Hooft interaction, and assuming that the coefficient of the anomaly term is temperature and 
density dependent. The special case of explicit breaking of chiral symmetry without Uy 4 (l) anomaly 
is also considered, since this scenario can provide additional information allowing to understand 
the interplay between the U^(l) anomaly and (spontaneous) chiral symmetry breaking effects. 
The pseudoscalar and scalar sectors are analyzed in detail bearing in mind the identification of 
chiral partners and the study of its convergence. We also concentrate on the behavior of the mixing 
angles that give us relevant information on the issue under discussion. In the region of temperatures 
(densities) studied, we do not observe signs indicating a full restoration of U(3)<8)U(3) symmetry 
as, for instance, the degeneracy of both oq and /o mesons with the pion. As we work in a real 
world scenario (m^ = << m^), we only observe the return to symmetries of the classical QCD 
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I. INTRODUCTION 


It is well known that Quantum Chromodynamics (QCD) has an approximate U(3)® U(3) 
chiral symmetry with its subsymmetry 11^(1) being explicitly broken by the axial anomaly 
In this context, the explicit and spontaneous breaking of chiral symmetry, as well as the 
Ua( 1) anomaly, play a special role, allowing for several nontrivial assumptions of low energy 
QCD: (i) the octet of the low-lying pseudoscalar mesons (vr, K, rj) consists of approximate 
Goldstone bosons; (ii) the rj — rj' phenomenology is characterized by large Okubo-Zweig- 
lizuka (OZI) violations. In fact, the important contribution of the U^(l)-breaking and the 
OZI rule violating terms, in the process of generation of meson masses and mixing angles. 


have been stressed in many phenomenological investigations 


New aspects of mixing 


and the consistent extraction of mixing parameters from experimental data have recently 


Q 


been discussed 

It is generally expected that ultra-relativistic heavy-ion experiments will provide the 
strong interaction conditions which will lead to new physics. In fact, it is believed that the 
availability of high-energy beams can provides the necessary conditions to observe small- 
distance scales, allowing to conhrm the QCD as the source of the strong interactions. 
Restoration of symmetries and deconhnement are expected to occur, allowing for the search 
of signatures of quark gluon plasma. 

The theoretical studies of QCD at hnite temperature and density present challenging 
questions, which may be the source of a productive complement for understanding relevant 
features of particle physics, not only in heavy-ion collisions, but also in the early universe 
and in neutron stars. In particular, the role played by the order of the chiral phase transition 
on the dynamical evolution of the systems, and possible experimental signs, have recently 


been addressed by some authors j^. In general, at hnite temperature and/or density one 
expects chiral symmetry to be restored above a certain temperature (density). 

In QCD, lattice calculations on the nature and order of the phase transitions indicate that 
light quarks experience a restoration of chiral ^mmetry as the temperature increases, with 
a transition temperature around 150 MeV 7|, |3,|9|. In the chiral limit, the restoration of 
chiral symmetry is signaled by the vanishing of the order parameters (qq) as the quark masses 
go to zero. The high temperature phase is sometimes described as a weakly interacting gas 
of quarks and gluons (plasma phase), which is clearly a simplistic picture for temperatures 
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around the transition temperature. It has been argued m that, just as in the more familiar 
low temperature phase, the behavior of the high temperature phase is characterized by the 
propagation of color-singlet objects. 

So far, the more reliable lattice QCD calculations for the phase transition have been 
focused on the non-zero temperature case. As an alternative to lattice QCD calculations, 
QCD-inspired models have been widely used in recent years to investigate hnite temperature 
and density effects. 

The assumption that the symmetric phase consists of mesonic modes and (deconhned) 
current quarks underlies the extended version of the Nambu-Jona-Lasinio (NJL) model 

nn 

111 I12| . This scenario allows to look for the spectrum of hadrons in parity doubling, whose 
degeneracy is taken as an indication of an effective restoration of chiral symmetries. In 
particular, scalar mesons and its opposite-parity partners, the pseudoscalars, are massive 
and degenerate in the symmetric phase. 

In the NJL model we can treat both the scalar and the pseudoscalar mesons on the 
same footing. The main problem concerning the scalar sector, = O’*", which has been 


under intense investigation over the past few years is that there are too many light 
scalars below 1 GeV. The two isoscalars /o(600) (a) an d fn f980) |l4| as well as the isovector 


ao(980) and the isospinor iLg(800) (that we will call k) [l^ scalars are enough candidates to 
£11 up a nonet of light scalars. Although it is accepted that large 4-quarks and meson-meson 


Q 


components jl^J are necessary to explain this nonet, here we shall assume a qq structure 
for the scalar mesons which are relevant to study the restoration of both chiral and axial 
symmetries. Recently, Dai and Wu [l^ claimed that (a, /o, oq, k) can be chiral partners of 
the pseudoscalar nonet (p,p',7r, K). Many other schemes have been suggested to describe 
the scalar meson properties. In fact, this is a very active field and no definitive conclusion has 
been reached as to which states are to be considered as qq, multi-quark, molecule, gluonia 


or hybrid states 


Q 


An important aspect of the problem is t 
symmetry in the restored chiral phase |l 


le role pla.v ed by the anomalously broken U^(l) 



JilLlU 


um, 


24l |. It has been argued that the 


chirally restored phase of QCD is effectively symme 


SU(Nj)(8)SU(Nj) at high temperature 


la |2ll 122, |23 


ric under U(Ny)®U(Nj) rather than 


25j |. Special attention has also been 


paid to whether or not the effective restoration of the Uyi(l) symmetry and the chiral phase 
transition occur simultaneously. This question is still controversial and is not settled yet. 
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indicating that we are still far from the full understanding of the dynamics of the processes 


under discussion. Here, we point out two scenarios discussed by Shuryak 


Q: 


m scenario 


1, Tc « 7[/(i) and the complete U(Nj)(8)U(Nj) chiral symmetry is restored well inside the 
quark-gluon plasma region; in scenario 2, ~ 

The effective restoration of the Uyi(l) symmetry means that all U^(l)-violating effects 
vanish, i.e., all order parameters of the U^(l) symmetry breaking must vanish. Since the 
origin of the anomalous interaction arises due to the presence of instantons in the physical 


state through the’t Hooft term 


26|, the effective restoration of the 11^(1) symmetry in the 


NJL model is equivalent to the vanishing of the effects of this interaction. 

The question is to look for observables which are strongly influenced by the anomaly and 
to see if they decrease and eventually vanish, indicating the absence of the anomaly. One 
of such quantities is the topological susceptibility, y, which, in purc color SU(3) theory, can 


be linked to the t]' mass trough the Witten-Veneziano formula j^|. The vanishing of this 
quantity could be an indication of the restoration of the Uyi(l) symmetry. In fact, lattice 
calculations at hnite temperature indicate a strong decrease of the topological susceptibility 


29|. and recent preliminary results at hnite density seems to conhrm this tendency 
30|. In addition, since the presence of the axial anomaly causes havor mixing, with the 
consequent violation of the OZI rule, both for scalar and pseudoscalar mesons, restoration 
of axial symmetry should have relevant consequences for the phenomenology of meson mixing 
angles, leading to the recovering of the ideal mixing. 


In a previous study m on effective restoration of chiral and axial symmetries in the 
NJL model, we have shown that the axial part of the symmetry is restored before the full 
U(3)®U(3) chiral symmetry. Here, we investigate two mechanisms to study an effective 
restoration of chiral and axial symmetries, which consists in two different ways for the be¬ 
havior of the coupling strength of the anomaly. One of them is based on a phenomenological 


ae coupiii 

HQ, 


decreasing nil 122 I, and the other one is inspired on the behavior of the topological suscep¬ 
tibility as indicated by lattice results at hnite temperature [^. This two cases are going to 
be compared with two limiting conditions: go = constant and Qd = 0 from the beginning. 
With this methodology we expect to disentangle the competition between U^(l) anomaly 
and chiral symmetry breaking effects. 

After the presentation of the model and the scenarios of restoration of the axial symmetry 
in Secs. H and HI, respectively, we start our investigation with the study of the consequences 
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of the effective restoration of chiral and axial symmetries with temperature and zero density 
(Sec. IV). Due to recent studies on lattice QCD at hnite chemical potential it is interesting to 
investigate also the restoration of the 11 ^( 1 ) symmetry at hnite density and zero temperature. 
In this case, we will consider two environment scenarios: completely symmetric matter 
{Pu = Pd = Ps) in Sec. V and quark matter simulating ’’neutron” matter in Sec. VI. Our 
conclusions are presented in Sec. VII. 


II. MODEL AND FORMALISM 


We consider the three-havor NJL type model containing scalar-pseudoscalar interactions 
and a determinantal term, the’t Hooft interaction generated by instantons in QCD, which 
breaks the Uyi(l) symmetry. The model has the following Lagrangian jhll ll^: 


C = - 7h)q+ ^gs ELo 

+ 9d {det [g (1 75 ) g ] -h det [g (1 - 75 ) g ] }• 

Here g = {u, d, s) is the quark held with three havors, Nf = 3, and three colors, W = 3. A“ 
are the Gell-Mann matrices, a = 0,1,..., 8, A*’ = -y/l I. 

Our ehective chiral held theory has the same chiral symmetry of QCD, coming out solely 
from quark interactions. The global chiral SU(3)®SU(3) symmetry of the underlying La¬ 
grangian (HD is explicitly broken by the current quark masses rh = diag(mu, m^, rris). As the 
Lagrangian ([Q) dehnes a non-renormalizable held theory, we introduce a cutoh which sets 
the 3-momentum scale in the theory. 


The NJL model can be generalized 
by applying the Matsubara technique 


:o the hnite temperature and chemical potential case 
33j | as is shown in the Appendix lA 21 


A. The gap equation 


In order to put the Lagrangian (1) in a form suitable for bosonization, it is useful to 
convert the six quark interaction in (HD into a four quark interaction [l^ |35|, ya |37| , 

allowing for the ehective quark Lagrangian: 


^eff = q{il^ - m)q 

+ 5'ab[(gA“g)(gA'’g)] -h Dab[(??75 A“g) (giqs A^)], 


( 2 ) 
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where the projectors Sab, Pab are presented in the Appendix lA II (Eqs. (jAlfl and ea). 

The bosonization procedure can be done by the integration over the quark helds in the 
functional integral with the effective Lagrangian (j2I), leading to an effective action (ESI) 
where, as shown in the Appendix lA 11 the natural degrees of freedom of low-energy QCD in 
the mesonic sector are achieved. 

The hrst variation of the effective action leads to the gap equation, 


Mi=mi- 2gg {qiqi) - 2g^ {q^qj) {qkqk) , ( 3 ) 

with i,j, k = u,d,s cyclic and Mi are the constituent quark masses. The quark condensates 
are determined by 

{mi) = -iM , = -iTr [^i(p)] , (4) 

where Si{p) is the quark Green function. 

B. Pseudoscalar and scalar meson nonets 

To calculate the meson mass spectrum, we expand the effective action (ESI) over the 
meson helds. Keeping the pseudoscalar mesons only, we hnd the meson masses by using the 
rest frame, P = 0, and the condition 

l-P,,n5(Po = M,P = 0) = 0. (5) 

For the non-diagonal mesons tt , A', we have 

P7r = 9s + 9d {qsqs) , ( 6 ) 


Pk = 9s + 9d {ququ) ■ ( 7 ) 

The polarization operator in Eq. (j^ takes the form given in the Appendix lA ll bv Eq. 


The quark-meson coupling and the meson decay constants /m are also evaluated according 


to the usual dehnitions 


Q. 


The inclusion of the’t Hooft interaction in the NJL model allows for havor mixing, giving 
rise to a P^-dependent mixing angle 6*p(P^) ^ ^ 0, |^. Our scheme for pseudoscalar 
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flavor mixing consists in the definition of the mixing angle 6p in such a way that 



where rj and rj' stand for the corresponding physical helds, and r]s and r]Q are the mathematical 
objects transforming as octet and singlet states of the SU(3)-flavor pseudoscalar meson 
nonet, respectively. 

The condition to diagonalize (Eq. (A4)) as = 

diag(P“^(P), gives us the equation for the mixing angle: 

2B 

tan26>p = ^ (9) 

as well as the inverse meson propagators, 


D-\P) = {A + C)- + 

d:Mp) = iA + c) + ^/{c-Ay + AB\ 


( 10 ) 

( 11 ) 


36 


38, 


39 


1^ . si 


since 


with ^ — Pgs — AIIoo(P), C — Poo — Aflgg {P),B — — (Pog + Allog(P)) and A — PooPgg — Pgg; 
the different projectors Pab and polarization operators nfj,(P) are dehned in the Appendix 
0(Eqs. (All-17)). 

In the rest frame, the condition P“^(Po = P = 0) = 0 and P“,^(Po = M^^f, P = 0) = 0 
gives as, as usual, the masses for the rj and r]'. 

As shown in other papers, in the framework of the NJL model 
A, B and C depend on P^, the mixing angles between the components rjo and rjs, Op (short 
notation of 6p{P‘^)) are P^-dependent. In the present paper, when studying temperature and 
density dependence of several quantities, we we only discuss the mixing angle for P^ = M^, 
for simplicity reasons; we checked that the behavior of the mixing angle for P^ = M^, gives 
information qualitively similar. 

The same technique used for the pseudoscalar sector can now be directly applied to the 
scalar resonances.We deal here with nine scalar resonances: three uq’s, which are the scalar 
partners of the pions, four k’s, being the scalar partners of the kaons, and the a and /o, 
which are associated similarly with the r] and r]'. As in the pseudoscalar case, we have mixing 
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between the a and /o and the neutral Oq . Keeping now the scalar mesons only, we have the 
effective meson action ()A2njl . The scalar meson masses are obtained from the condition 


i-s'«nf,(p„ = Af,p = o) = o, 


( 12 ) 


with 


9s 9D {QsQs) ) 


(13) 


Sk = 9s- 9d {ququ) ■ (14) 

The polarization operator is presented in the Appendix I A 11 Eq. (IA21|1 . 

Finally we can determine the meson masses of oq and k using the respective dispersion 
relations: 


s'„„ni(M,„,o) = 0, 

(15) 

■^.nf,(M.,o) = o. 

(16) 


For the diagonal mesons Og, a and /g we take into account the matrix structure of the 
propagator in (jA20j) . In the basis of Og — a — /o system, we write the projector Sab and the 
polarization operator Ilfj, as matrices (see the Appendix To find the masses of the a 

and /o mesons we use the inverse propagator of the corresponding mesons as indicated in 
the Appendix lA 11 The value of the angle 6s can also be fixed by a condition similar to Eq. 

When Pg > Mi + Mj, i. e., when the mass of the meson exceeds the sum of the masses of 
its constituent quarks, the meson can decay in its quark-antiquark pairs, being, therefore, a 
a resonant state. Then, Eqs. ®, m and m have to be calculated in their complex form 
in order to determine the mass of the resonance Mm and the respective decay width Tm- 
Thus, we assume that these set of equations has solutions of the form 

Pg = Mm — (17) 

and, on the other hand, we have to take into account the imaginary part of the integrals 
(A8) (for details see Appendix Al). 
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C. Vacuum properties and model parameters 


The NJL model exhibits a vacuum phase where chiral symmetry is spontaneously broken, 
a mechanism which generates the constituent quark masses. The model is fixed by the 
coupling constants gs, Qd the Lagrangian (^Q), the cutoff parameter A which regularizes 
momentum space integrals !{ and and the current quark masses rrii. We start by 

considering two sources of chiral U(3)®U(3) symmetry-breaking: (i) current quark masses; 
and (ii) Uyi(l) symmetry-breaking effective interaction. 

As already referred, the SU(3) version of the NJL model has five parameters, and we would 
expect a priori that one can uniquely fix those parameters in order to fit five observables f.^ , 
AAr, Mx, and M^j'. However, this is not the case as can be seen comparing the parameter 
sets of [11] and [12]. We follow the methodology of Ref. [12] and set rriu to the value 5.5 
MeV, and fix the remaining four parameters by fitting /^r, AAr, Mk and Mr,'. The rj meson 
in this way is predicted with a mass of 514.8 MeV. This allows for a good overall agreement 
of our numerical results with the experimental or phenomenological quantities as shown in 
Table I. 

However, we point out that this prescription has some problems in which concerns to the 
description of the rj' meson. As is well known the NJL model does not confine. Formally, 
this is reflected by the fact that integrals like and hence the polarization function for 

some mesons, get an imaginary part above the gg-threshold that is calculated as indicated 
at the end of the previous section. 

We will consider a second parametrization without U^(l) symmetry-breaking effective 
interaction {qd = ^) which is also presented in Table 1. With this parametrization we also 
have an overall satisfactory fit to meson properties and quark condensates at zero tempera¬ 
ture and density. However, as expected, the results show that the anomaly term is necessary 
to obtain the correct meson mass spectra, especially by giving the rj' and Oq its large masses, 
as well as the splitting between n/r], and a/ao meson masses. 
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Explicit symmetry breaking with U/i(l) anomaly (qd ^ 0) 


Parameter set 

Physical quantities 

and constituent quark masses 

U = 92.4 MeV 

niu = rud = 5.5 MeV 

= 135.0 MeV 

rus = 140.7 MeV 

Mk = 497.7 MeV 

A = 602.3 MeV 

M^/ = 960.8 MeV 

gsA^ = 3.67 

= 514.8 MeV* 

goA^ = -12.36 

fx = 97.7 MeV* 

Mu=Md = 367.7 MeV* 

= 728.8 MeV* 

Ms = 549.5 MeV* 

Mao = 873.3 MeV* 


= 1045.4 MeV* 


M/o = 1194.3 MeV* 


0p = -5.8°* ; 03 = 16°* 


Explicit symmetry breaking without Ua( 1) anomaly {go = 0) 

U = 92.4 MeV 

niu = md = 5.5 MeV 

= Mn = 135.0 MeV 

nis = 138.75 MeV 

Mk = 497.7 MeV 

A = 602.3 MeV 

fx = 95.4 MeV* 

gsA^ = 4.64 

Mr,' = 707.5 MeV* 

goA^ = 0 

M^ = Mao = 740.1 MeV* 

Mu= Md = 368 MeV* 

M^ = 985.38 MeV* 

Ms = 587.4 MeV* 

M/o = 1194.8 MeV* 


Op = -54.74°* ; 0s = 35.264°* 



TABLE I: Physical quantities in the vacuum state and the parameter sets for the two symmetry 
breaking patterns studied in this work. The asterisk signalize predicted physical quantities. 
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III. SCENARIOS OF RESTORATION OF THE AXIAL SYMMETRY AND EN¬ 
VIRONMENT CONDITIONS 


Model calculations, for instance within NJL type or sigma models, generally describe 
the restoration of chiral symmetry as a natural consequence of the increase of temperature 

es associated to the anomaly, although 


or density. However, it is found that the observab 

nn 

decreasing, do not show a tendency to vanish [SlJ, |^, 


42 


43| |. The anomaly in our model 


is present via the’t Hooft interaction and its effects appear explicitly in the gap equations 
© and in the mesons propagators through products of the anomaly coefficient by quark 
condensates (see the expressions of the projectors Sab, Pat given by Eqs. (EH) and (El)- 
Such quantities, that act as a kind of’’effective anomaly coupling”, will be denoted from now 
on as {gD)i = go {gi^i)- The vanishing of such effective coupling should imply the vanishing 
of the observables associated to the anomaly. However this does not happen in the present 
model without being enforced because, while the non-strange quark condensates decrease 
asymptotically, leading to an effect almost negligible of ((s'n))^)) fh® same does not 

happen with {go)gi since restoration of chiral symmetry does not occur in the strange sector 
and {qsqs) has always an appreciable value. Therefore, the vanishing of {gD)i, in general, 
should be accomplished by assuming that the anomaly coefficient is a decreasing function 
of temperature or density. 

More attention has been paid, up to now, to the restoration of axial symmetry with 


Dy the lattice results for 
, that indicate a consid- 


temperature than with density, a motivation which is supported 
the behavior of the topological susceptibility with temperature j2^ 
erable decrease of this quantity. However, the theoretical arguments concerning the possible 
restoration of axial symmetry, whether temperature or density are considered, are similar. 
While the difficulties in testing the QCD vacuum at high density in heavy-ion collisions 
are not yet removed, but expecting that this will hopefully happen in future experiments, 
it is useful to have predictions for the non-perturbative regime, even at a qualitative level. 
Model calculations in NJL model, although not being an alternative to lattice calculations, 
can provide a useful contribution. Moreover, lattice calculations for the behavior of the 


topological susceptibility with density [^, although still in a early stage, suggest that this 
observable is also a decreasing function of density. In view of the considerable interest in 
the investigation of the behavior of matter at high densities, and the possible restoration of 
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Anomaly coefficient 

Case I 

Case II 

Case III 

Constant 

Fermi function 

Decreasing exponential 


TABLE II: Different schemes of explicit axial symmetry breaking with temperature (density). 


symmetries under these conditions, it is certainly worthwhile to do an exploratory study on 
the restoration of the axial symmetry by assuming that go is density dependent, in a form 
similar to the temperature dependence. 

So, after considering the extreme case of a constant anomaly coupling, go, we will consider 
2 scenarios to study the effective restoration of axial symmetry as summarized in Table II. 

Case I: the anomaly coefficient go constant for all range of temperatures or densities. 

Case II: the anomaly coefficien t Qp is a dropping function of temperature or density. 
Following the methodology of Ref. jd^, the temperature dependence of go is extracted by 
making use of the lattice results for the topological susceptibility, y, . The expression for 
X in the NJL model is presented in the Appendix IA 81 Eq. (jA42jl . In view of the arguments 
presented above, it seems reasonable to model the density dependence of Qn extrapolating 
from the results for the finite temperature case and proceeding by analogy j31|. 

Case III: the anomaly coefficient has the form of a decreasing exponential {goiT) = 
5'£)(0)exp[—(T/T q)^]). This phenomenological pattern of restoration of the axial symmetry 


was proposed by Kunihiro ^ in the framework of the present model. Here we consider 
a dependence of the anomalous coupling constant on density also inspired on the finite 
temperature scenario. 

We also consider a simplistic scenario without U^(l) anomaly (go = 0), which is achieved 
in our model by choosing the second parametrization presented in Table I. We expect that 
this scenario, being a limiting case, might provide additional information allowing to un¬ 
derstand the interplay between the U^(l) anomaly and flavor symmetry breaking effects. 
In fact, in this case the dominant effects come from spontaneous chiral symmetry breaking 
through quark loop dynamics. 

For a more complete understanding of the density effects we will consider two differ¬ 
ent scenarios of quark matter: (i) symmetric quark matter; and (ii) neutron matter in 
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FIG. 1: Temperature dependence of {gD)u panel) and {go)s (right panel) for the different 


cases. 


/9-equilibrium. So, the different patterns of axial symmetry with Qd = constant (Case I), 
Case II and Case III, and go = 0, are going to be studied in hot media, in symmetric quark 
matter and in neutron matter. 

The restoration of chiral symmetry with temperature or density has been extensively 


studied in the present model with go constant [Ij, |38|, 



45| |. A general conclusion of such 


studies is that chiral symmetry is effectively restored in the SU(2) sector, but, in the range 
of densities or temperatures generally considered, the same does not happen in the strange 
sector. It should be noticed that, as we will show, this conclusion will not be affected by 
the different patterns of axial symmetry restoration here considered. 

Since in all cases chiral symmetry is explicitly broken by the presence of non-zero cur¬ 
rent quark mass terms, chiral symmetry is realized through parity doubling rather than by 
massless quarks. So, the identihcation of chiral partners and the study of its convergence is 
the criterion to study the effective restoration of chiral and axial symmetries. 


13 
























IV. RESULTS FOR THE MESONIC BEHAVIOR AT FINITE TEMPERATURE 
AND ZERO DENSITY 

In this section we analyze the mesonic behavior at hnite temperature and zero chemical 
potentials. A signihcant feature of this analysis is that the observables, that depend on 
the anomaly coupling only via {gD)u ({ 90 ) 0 )^ signihcantly affected by the specihc 

temperature dependence of go, in the high temperature region, because chiral symmetry is 
approximately restored with the consequent asymptotic vanishing of the non-strange quark 
condensates. In order to see the importance of the behavior of the effective anomaly coupling 
for the quantities under study, we plot them in Fig. [TJ 


A. Explicit chiral symmetry breaking with Uy4(l) anomaly 


Case I. We will start with Case I (see Fig. a) and (b)) that will be compared with 
the other cases under discussion. In the panel (b) we have the meson masses for the K- 
meson and its chiral partner k. In the panel (a) we plot the other scalar and pseudoscalar 
mesons: (cr, /o, Oq, tt, rj, rj'). In both panels, and for all graphics, the dotted line means the 
respective continuum. The crossing of the tt and rj lines with the quark threshold 2M„, and 
the K line with Mu + Mg indicates the respective Mott transition temperature, Tm- Mott 
transition comes from the fact that mesons are not elementary objects but are composed 
states of qq excitations, and is dehned by the transition from a bound state to a resonance in 
the continuum of unbound states. Above the Mott temperature we have taken into account 
the imaginary parts of the integrals and used a hnite width approximation 

Let us summarize here the behavior of the pseudoscalar mesons and analyze what this 
behavior can tell us about possible restoration of symmetries. One can see that Mott tem¬ 
peratures for rj and vr mesons are: Tm^ = 180 MeV and Tm^ = 212 MeV. The vr and K 
mesons become unbound at approximately the same temperature: Tmj^ =210 MeV. On the 
other side, the g' is always above the continuum uju = 2Mu, and g has always a strange 
component for all temperatures, once its mixing angle 6p never gets the ideal value: 6p = 
-54.736° (see Fig. E)). 

Concerning the scalar sector, we notice that the a meson is the only scalar meson that 
is a bound state for small temperatures (the others are always resonances) but turns into 
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FIG. 2: Temperature dependence of meson masses and of the continuum thresholds (dotted lines) 
2 Mu, 2Ms (left panels) and Mu + Mg. 


a resonance at Tm^ — 160 MeV. This meson has a strange component that decreases with 
temperatnre bnt never vanishes since the ideal mixing angle, 6s = 35.264°, is never attained 
in the range of temperatnres stndied (see Fig. El) 

For T > 250 MeV the a starts to be degenerate with the vr. As for the oq meson, it is 
always a non-strange state and is above the continnnm Uu = 2Mu. It can be seen in Fig. 
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FIG. 3; Left panel: scalar and pseudoscalar mixing angles as a function of temperature for the 
three cases presented in Table II, and for the ideal mixing. Right panel: Topological susceptibility 
as a function of temperature for the three cases. The lattice data results are plotted with error 
bars 
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n a) that the partners (vr, a) and (?], oq) become degenerate at almost the same temperature. 
In both cases, this behavior is the signal of the effective restoration of chiral symmetry in the 
non-strange sector. Distinctly, the rj' and /o masses do not show a tendency to converge in 
the region of temperatures studied. We interpret this behavior as an indication that chiral 
symmetry does not show tendency to get restored in the strange sector (see 2 Ms, upper 
dotted curve in Fig. l2Ka)). 

Finally, we focus on the k- meson (Fig. Etb)). It is always an unbound state and, as the 
temperature increases, it shows tendency to get degenerate in mass with the iF-meson. For 
comparison purposes, we summarize in Table III an overview of the transition temperatures 
of the effective restoration of chiral (second line) and axial (third line) symmetries in the 
different cases studied in the present paper. The masses of the corresponding chiral partners 
become degenerate above the referred temperatures. 

Summarizing, the SU(2) chiral partners (tt, a) and {r],ao) become degenerate at T ~ 250 
MeV; the chiral partner {K,k) converges at T ~ 350 MeV and [r]', /o) do not show a 
tendency to converge in the region of temperatures studied. 

As expected, the axial symmetry is not restored at high temperatures and the topological 
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Case I 

idD = constant) 

Case II 

Case III 

o 

II 

Q 

SU(2) chiral- 

transition temperature 

250 MeV 

250 MeV 

225 MeV 

300 MeV 

U(2) axial- 

transition temperature 

— 

350 MeV 

225 MeV 

— 


TABLE III: Transition temperatures of the effective restoration of chiral and axial symmetries in 
the different cases. 


susceptibility is also far away from being zero (see Fig. ED- 

Case II. Some of the results for Case II have been presented in j31|. Here we summarize 
the conclusions obtained. 

As rriu = rrid « the (sub)group SU(2)®SU(2) is a much better symmetry of the 
Lagrangian (1) than SU(3)0SU(3). So, the effective restoration of the SU(2) symmetry 
implies the degeneracy between the chiral partners (vr, a) and (p, oq) which is verified around 
T ~ 250 MeV (see Fig. Efc) and Table III). For temperatures at T ~ 350 MeV both oq 
and a mesons become degenerate with the tt and r] mesons, showing, as explained below, an 
effective restoration of both chiral and axial symmetries. Without the restoration of U^(l) 
symmetry (Case I), the oq mass was moved upwards and never met the vr mass, the same 
argument being valid for the a and rj mesons. We remember that the determinant term acts 
in an opposite way for the scalar and pseudoscalar mesons as can be seen, for instance, in 
Eqs. (0 and (USD. So, only after the effective restoration of the U^(l) symmetry we can 
recover the SU(3) chiral partners (tt, an) and (? 7 , n) which are now all degenerate. This is 


compatible with scenario 1 of Shuryak 1^: the signals for the effective restoration of the 
axial symmetry occur at a temperature where the signals of the full restoration of U(3)(8)U(3) 
symmetry are not yet visible. In fact, the rj' and /o masses do not show a clear tendency to 
converge in the region of temperatures studied, this absence of convergence being probably 
due to the fact that, in the region of temperatures above T ~ 350 MeV, those mesons are 
purely strange and the chiral symmetry in the strange sector is far from being effectively 
restored. 

The analysis of the temperature dependence of the mixing angles in Fig. El allowing 
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for a better understanding of the meson behavior through the evolution of the quarkonia 
content, provides further indication of the restoration of the axial symmetry: 6s {Op) starts 
at 16° (—5.8°) and goes, smoothly, to the ideal mixing angle 35.264° (—54.74°). This means 
that flavor mixing no more exists. In fact, referring to the SU(2) chiral partners (vr, a) and 
{ri,ao), we found that the Oq and vr mesons are always purely non-strange quark systems, 
while the a [r]) meson becomes purely non-strange when 6s {6p) goes to 35.264° (—54.74°), 
at T ~ 350 MeV. 

Analyzing the {K,k) chiral partner, we conclude that the behavior of the mesons is not 
signihcantly influenced by the type of temperature dependence of gp used here, as expected. 
In fact, in the range of temperatures where goiT) could be important. Mg does not change 
appreciably, and we know that these meson masses are very sensitive to Mg. For the range 
of temperatures where the (tt, uq) and {rj, a) chiral partners become degenerate, the strange 
quark mass of Mg is already independent of the gp dependence of the temperature (see 
2Ms, upper dotted curve in Fig. El^c)). This is due to the fact that, as explained before. Mg 
depends on the anomaly through {gp)^- 

We notice that our analysis of the effective restoration of symmetries is based on the 
degeneracy of chiral partners that occurs in a region of temperatures where the mesons are 
no more bound states (they dissociate in qq pairs at their respective Mott temperatures 
ji^ l^h Moreover, the mesons g' and /o are qq resonances from the beginning and its 
description is unsatisfactory. 

Summarizing, we conclude that at T ~ 250 MeV the SU(2) chiral partners become 
degenerate in mass, whereas at T ~ 350 MeV, the same happens with (tt, a, g, oq) mesons: 
the OZI rule is restored and y goes asymptotically to zero (Fig. El dashed line of right 
panel). These results indicate an effective restoration of the Ua( 1) symmetry. 

Case III. Finally, we analyze Case III that is similar to Case II as we can see in Fig. 
121 The main difference is that the temperature dependence of gp used does strengthen 
signihcantly the chiral phase transition. In fact, the SU(2) chiral partners (tt, a) and {g,ao) 
are all degenerate for T ~ 225 MeV (T ~ 250 MeV in Case II). Linking this fact to the 
behavior of the y (in Fig. El dotted line in right panel) that goes very fast to zero, being 
zero at about 250 MeV, and with the behavior of the mixing angles (Fig. El dotted lines in 
left panel), 6 p and 6 s, that go both to its ideal values at 200 MeV, we conclude that both 
symmetry restorations happen around the same temperature. 
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The comparison between Case \ {qd = constant) and Case III helps to understand this 
situation. We observe that the more rapid decrease of the temperature dependence of Mg in 
Case III, till T ^ 250 MeV, indicated by the upper dotted lines {2Ms) in the left panel of Fig. 
IKsee Fig. |2Da), (e)), cooperates with the decreasing of goiT) allowing for the restoration of 
chiral and axial symmetries at the same temperature T ~ 225 MeV. This can also be seen, 
for instance, in the behavior of the effective anomaly coupling {gD)s goes to zero at 
almost the same temperature (see Fig. Q). This is in accordance with scenario 2 of Shuryak 


Q. 


The existence of cooperative effects of restoration of chiral and axial symmetries has 


already been noticed by Kunihiro who report a situation where the axial symmetry is 
restored before chiral symmetry, a scenario usually considered not realistic . 

Concerning {K,k) chiral partner, the conclusions are similar to those of Case II; the only 
difference is a faster decrease of the splitting in the low temperature region, due to the faster 


decrease of Mg. 


B. Explicit chiral symmetry breaking without Ua( 1) anomaly 

We consider now a simplistic scenario without Ua( 1) anomaly (go = 0), which is achieved 
in our model by choosing the second parametrization presented in Table I. 

We start with the vr-meson that, as expected, is always degenerate with g. In fact, the 
meson is a pure non-strange state for all temperatures, with a ideal mixing angle 9p = 
-54.736°. 

On the other side, the g' is always a pure strange state which crosses the continuum 
uju = 2Mu for T > 110 MeV, becoming then a resonance state. Like in the gp ^ 0 cases, the 
g' meson shows no tendency to become degenerate with /q, a consequence of the insufficient 
restoration of chiral symmetry in the strange sector, as it has already been noticed. The 
ao (= cr) is always a non-strange state, it is always above the continuum oju = 2M„. 

As the temperature increases, due to the absence of the 11^(1) anomaly, the members 
of the chiral pairs (tt, a) and {g,ao) become all degenerate simultaneously (T ~ 300 MeV), 
reflecting the effective restoration of chiral symmetry in the non-strange sector. We no¬ 
tice that, as indicated in Table III, this is the case where the transition temperature to 
the SU(2)®SU(2) symmetry is higher, indicating that, as already referred, the anomalous 
coupling can be important to drive the effective restoration of the chiral symmetry itself. 
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Concerning the kaon and its chiral partner k, they show a clear tendency to get degener¬ 
ate, but at temperatures that are higher than in the previous cases. 

Summarizing, the high temperature regime (T 300 MeV) in Cases II and III, where 
the axial symmetry is effectively restored, and the situation Qd = ^ are very similar: the 
SU(3) chiral partners (tt, oq) and (rjjCr) are degenerate, and the rj' and /o mesons have 
similar splittings. The more relevant differences in the behavior at lower temperatures are 
manifestations of the different role played by the axial anomaly and the dynamical flavor 
symmetry breaking effects. For instance, the constituent strange quark mass has a very 
different behavior in the three scenarios as can be seen by the curve representing 2 Mg 
(upper dotted curve) in Fig. |21 ((c), (e) and (g)). In addition, we notice that, differently 
from Case II, in Case III the restoration the U^(l) symmetry drives of chiral symmetry. 


V. RESULTS FOR SYMMETRIC QUARK MATTER 


In order to study the effective restoration of chiral and axial symmetries at finite density, 
we start by considering a completely symmetric quark matter (p„ = Pd = Ps)- Before 
our analysis, let us make some considerations about this type of matter. Although rather 
schematic, this case simulates a situation where the hypothesis of absolutely stable strange 
quark matter (SQM) can be explored j36|. It has been argued that SQM may only 

be stable if it has a large fraction of strange quarks {ps ~ Pu ^ Pd)- The speculations on 
the stability of SQM are supported by the observation that the inclusion of the strangeness 


degree of freedom allows for a 
a sizable binding energy. In 


B er decrease of the strange quark mass which can produce 
we have confirmed this tendency when compared with 
neutron matter. We notice that there are always strange valence quarks present, so the 
strange quark mass decreases more strongly, although, even in this case, it is still away 
from the strange current quark mass for high densities a. TKe advantage of—, 
this type of matter is that, like in the non-zero temperature case, all three pions and all 
four kaons are degenerate in medium (contrarily to what happens for neutron matter in 
/d-equilibrium to be discussed in the next section). So, the present environment can provide 
fruitful comparisons with the non-zero temperature case. 
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P /Pn P /Pn 

' sym ' 0 ' sym ' 0 


FIG. 4: Density dependence of meson masses and of limits of the Dirac sea continua (dotted lines) 
defining qq thresholds for the mesons r]', oq, k. The /o meson is always a resonance state. 

A. Explicit chiral symmetry breaking with Uyi(l) anomaly 

The study of density effects is performed using a methodology analogous to the temper¬ 
ature case. So we will also consider 3 scenarios for the behavior of the anomalous coupling 
qd (see Table II). Some conclusions are very similar to the temperature case as we can check 
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in Fig. EJ so we will concentrate on the main differences. 

Case I. In this case we have an interesting phenomenon: the oq does not degenerate 
with the 77 -meson, but with the r^'-meson as we can see in Fig. EJa). In fact, the r^'-meson, 
that starts as an unbound state and becomes bound for psym > 4.5po, degenerates with the 
oo-meson for higher densities. This is due to the presence of strange valence quarks in the 


36j |. so the influence of the s sector 


medium, which causes Mg to decrease more strongly 
will be lower in the mass of the r]'. 

Case II. The study of Case II in symmetric quark matter is inspired, as already referred, 
in the previous Case II at finite temperature. So, we postulate a dependence for y formally 
similar to the temperature case as is shown in Fig. El (dashed line), i.e., using a Fermi 
function. With this topological susceptibility we obtain the density dependent anomalous 
coupling gniPsym)- 

Using this density dependence we arrive at conclusions similar to the finite temperature 
case discussed before. The chiral partners (tt , a) become degenerate at a density psym — 
3.5po and the same happens to the chiral partners ( 77 , 00 ) (Fig. life)); this density is the 
onset for effective restoration of chiral symmetry in the SU(2) sector. 

The analysis of the mixing angles (Fig. Ej) indicate that at psym — 4po the scalar and 
pseudoscalar mixing angles reach its ideal values and, consequently, the 77 and the a become 
purely non strange. At this density the 77 ' becomes purely non strange and does not show a 
tendency to degenerate with /o, as in the finite temperature case. 

Summarizing, as the density increases, the chiral partners (vr, a) and ( 77 , oq) become de¬ 
generate (for Psym — 4po)- Associating this with the behavior of the 6 p and 63 mixing angles 
and the behavior of the chiral susceptibility, that goes to zero (dashed line in Fig. El), we 
conclude that we have an effective restoration of Ua( 1) symmetry in this situation. 

Case III. Similarly, for Case III we postulate the following dependence for gp : 
gniPsym) = 5 'D( 0 )exp[—(psym/po)^], wliich is inspired in the corresponding finite tempera¬ 
ture scenario. The topological susceptibility with this dependence of the coupling anomaly 
is plotted in Fig. El dotted line. This case is very similar to Case II and the overall conclu¬ 
sions are parallel to the finite temperature case (to compare see Fig. Efe), (f)). The density 
dependence for gniPsym) that we used also does strengthen the chiral phase transition: both 
symmetry restorations happen simultaneously for slightly lower densities {psym — S.Opo)- 

In fact, the SU(3) chiral partners (tt, oq) and ( 77 , a) are all degenerate at very earlier values 
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FIG. 5: Left panel: scalar and psendoscalar mixing angles as a fnnction of density for the three cases 
presented in Table II, and for the ideal mixing {g£) = 0). Right panel: topological susceptibility as 
a function of density for the three cases. 

of the density [psym, — S.Opo), compared with Case II {psym — 4.0po)- This results from the 
behavior of x (in Fig. El dotted line in right panel) that goes to zero for psym — ^Po, and by 
the behavior of the mixing angles (Fig. El dotted line in left panel) where both, 6 p and 6 s, 
go to the ideal mixing angles for psym — 2.5po- 

Concerning the {K,k) partners (Fig. Eff)), we conclude that in all three cases their 
behavior is very similar: they practically do not depend on the shape of gp. 

B. Explicit chiral symmetry breaking without Ua( 1) anomaly 

Finally, we analyze the behavior of mesonic modes in the absence of the axial anomaly: 
gD = 0. Looking at the 7r-meson behavior plotted in Fig. 111(g), we conclude that the pion 
is always degenerate with rj and they are always bound states. The rj {rj') meson is a pure 
non-strange (strange) state for all range of densities. For 2po ^ Psym ^ 4 po fhe p'-meson is 
a resonance state as can be seen in Fig. Hg)- As in the Cases I, II and III, the p'-meson 
does not show tendency to become degenerate in mass with the /o-meson. 

The Oo (cr) is always a non-strange state and, for psym = 0, its mass is higher than 
uju = 2M„. As the density increases, it immediately becomes a bound state and degenerates 
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with its chiral partner rj (vr) for densities psym — 5po- So, for psym ^ 5po the four bound state 
mesons (vr, a, p, Oo) become degenerate reflecting the effective restoration of chiral symmetry 
in the non-strange sector. 

Concerning the kaon and its chiral partner k (Fig. 0(h)) they show a clear tendency to 
get degenerate at high densities, where both mesons are bound states. We remark that the 
degeneracy of chiral partners in symmetric quark matter occurs in regions where the mesons 
are bound states. 

Finally, the more signihcant difference between case Qd = ^ and the other cases is that 
the chiral symmetry effective restoration occurs latter, similarly to the situation at non-zero 
temperature. 


VI. RESULTS FOR ’’NEUTRON” MATTER IN (3 -EQUILIBRIUM: 


We consider now asymmetric quark matter in weak equilibrium and charge neutrality, 
supposedly of the same type of the existing in the interior of neutron stars. To insure this 
situation, we impose the following constraints on the chemical potentials and densities of 
quarks and electrons: 

^^d h's T h'e 


and 


2 

gPu 


■^{Pd + Ps) — Pe 


0 , 


(19) 


with 

Pi = ~ - ^i) and Pe = /ie/37r^ (20) 

Similarly to the hnite temperature case, and as already explained in the section III, 
chiral symmetry is effectively restored only in the SU(2) sector, in the range of densities 
considered, a conclusion that is independent of the specific form of the dependence on 
density of the anomaly coefficient, gr,. The effective anomaly coupling, shown in Fig. 0 
although exhibiting details different from the hnite temperature and from the symmetric 
quark matter cases, are qualitatively similar. 

Let us emphasize some specihc aspects on the behavior of the strange quark mass with 
density. Although in the present case, at low densities, there are no strange quarks in the 
medium, the mass of the strange quark decreases, although smoothly, due to the effect of 
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FIG. 6: Density dependence of {90)^ (1®^^ panel) and {go) g (I'ight panel) for the different cases. 

the’t Hooft interaction; eventnally it becomes lower than the chemical potential for strange 
qnarks (at ps — 3.8po)- A more prononnced decrease of the strange qnark mass is then 
observed, which is no more due to the anomaly (we can see from Fig. IHlthat (s'd)^ is already 
very small) but to the presence of valence strange quarks in the medium (see Eq. fj2()j) h 

Concerning the meson spectra and the mixing angles, we will show that new aspects also 
appear, mainly in the high density region, and will be discussed in the sequel. 

As it is the only scenario where the flavor symmetry exhibited by the physical vacuum 
state is violated by the weak interaction conditions (HHD, this implies several consequences: 

(i) splitting between charge multiplets of pions and kaons; 

(ii) appearance of low-lying modes above a certain density. 

This leads us to focus on the behavior of all nine pseudoscalar mesons and respective 
scalar partners, as well as on the chiral partners of the low-lying excitations. Before we start 
our discussion we remark the following: 

1. We start by analyzing the chiral asymmetry parameter which is a measure of the 
violation of the isospin symmetry. 

2. We will follow the structure used for the hnite temperature and completely symmetric 
matter which leads to the study the scenarios: Cases I, II and III, with qd 7 ^ 0; and 
the case go = 0 . 
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FIG. 7; Chiral asymmetry parameter as a function of density for the three cases presented in Table 
II and the case qd = 0 - 


A. Chiral asymmetry parameter 


Solving the gap equation © one verifies that, in the different cases summarized in Table 
II, the constituent quark mass decreases slightly more than as the density increases. 
Bearing in mind a qualitative analysis of the effects of chiral symmetry breaking (restoration) 
it is useful to plot the isospin asymmetry parameter 


\M^-Md\ 
+ M, ’ 


( 21 ) 


as a function of the baryonic density in the several cases under discussion. As is shown in Fig. 
[Tithe chiral asymmetry parameter xa is more signihcant in the absence of the anomalous 
coupling constant. We remark that the presence of the anomaly in the model has the effect 
of reducing the isospin asymmetry in an SU(2) broken system like the neutron matter case. 
The main consequences of this isospin asymmetry of the medium must be visible in the 
behavior of chiral partners. 


B. Explicit chiral symmetry breaking with Uyi(l) anomaly 


Case I. T he p seudoscalar sector in neutron matter in /^-equilibrium was extensively 


studied in |34, 


35 


3 , 


which corresponds to Case I of the present work. Here we will focus 
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Pb^Po Pb^Po 


FIG. 8: Density dependence of meson masses and of limits of the Dirac sea continua (dotted lines) 
defining qq thresholds for the mesons. The low-lying solutions are also included. The anomaly 
coupling constant is kept constant (Case I). 

mainly on the possible degeneracy of chiral partners. 

In Fig. a), the meson masses are plotted as functions of the density. The SU(2) chiral 
partners (tt^, a) are always bound states. The pion is a light quark system for all range 
of densities and the a meson has a strange component at = 0 but never becomes a 
purely non-strange state because 9s never reaches 35.264°, the ideal mixing angle (Fig. IHl 
Case I). As the density increases these mesons become degenerate (ps > 3po). At the same 
density, the SU(2) chiral partner (p, oq) is also degenerate. The p-meson is always a bound 
state, contrarily to Oq that starts as a resonance, once its mass is above the continuum, and 
becomes a bound state for ps > 0.5po. However, the Oq mass separates from the p mass 
and goes to degenerate with the p'. To understand this behavior we need to look for the 
behavior of the mixing angle 6p. From Fig. EJ Case I, we observe that the angle 6*p, which 
starts at —5.8°, changes sign at pp ~ 3.5po becoming positive and increasing rapidly, which, 
as will be seen, we interpret as an indication of a change of identity between p and p'. 
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We remember that up to the density ps = 3.8po (see end of subsection IVI Afl the {qq)s = 
ss content is induced by the mixing effects only. Above this density strange valence quarks 
are present (see Eq. PH)) and induce the strange quark mass to decrease faster. 

This behavior induces changes in the percentage of strange, {qq)s = ss, and non-strange, 
= \{uu + dd), quark content in rj and p' mesons: at low density, the p' is more strange 


than the p, but the opposite occurs at high density 


35j |. Then p' will degenerate in mass 


with the oo-meson that is always a non-strange state. Finally, the /o-resonance is always a 
strange state that shows no tendency to become degenerate with any other meson. 

Now let us comment on the tt^ behavior and the respective chiral partners that are 
plotted in Fig. HKb). The vr^-mesons are always bound states and their masses increase with 
density. On the other side, the Of^-mesons start as resonances and become bound states: 
the Oq” at pB — 0.25po, and the at ps — O.Spo- However, they never degenerate with the 
respective pions in the considered range of densities. This is, once again, due to the fact that 
the chiral symmetry in the strange sector is not restored, and the absence of the mechanism 
of restoration of the Ua( 1) symmetry is also relevant in this context. This will influence 
the behavior of the tt^ and mesons through Eqs. © and (HI, respectively, because the 
quark condensate (ss) is still very high (see Fig. El for (s'd)^)- A different scenario occurs 
for kaons and their chiral partners: and in Fig. IHtc) and K^{K^) and kP{kP) in Fig. 

Hd). 


Q, 


As is has already been shown Ea, below the lower limit of the Fermi sea continuum 
of particle-hole excitations, there are low bound states with quantum numbers of K~ , 
and vr’*'. Here we show that these low-energy modes, collective particle-hole excitations of 
the Fermi sea, have corresponding chiral partners. The behavior of the low-energy chiral 
partners with density is similar to that of the respective high energy modes and does not 
present meaningful differences in Cases I, H and HI. This can be seen in Fig. H c) and (d), 
for {K~, K~) and (A'°,k°), respectively. 


We also saw in m that in the present approach the criterion for the occurrence of kaon 
condensation is not satished since the antikaon masses are always larger than the difference 
between the chemical potential of strange and non-strange quarks. This conclusion is still 
valid in the other cases. 


Case II. Like in the previous section, we postulate in Case H a density dependence of y, 
as a Fermi function, formally similar to the hnite temperature case (see Fig. IHl dashed line 







FIG. 9: Left panel: scalar and pseudoscalar mixing angles as a function of density for the three 
cases and for the ideal mixing. Right panel: topological susceptibility as a function of density for 
the three cases. 

in the left panel). Then we can apply this dependence to model the anomalous coupling, 
allowing the calculation of all observables. 

Analyzing the mixing angles (Fig. IHl right panel) we observe that the behavior of 6s is 
similar to the non-zero temperature and completely symmetric quark matter cases: it starts 
at 16° and increases up to the ideal mixing angle 35.264°. A different behavior is found for 
the angle 9p, that changes sign at pp — 4.9po (— 3.50po in Case I): it starts at —5.8° and 
goes to the ideal mixing angle 35.264°, which also leads, by similar reasons as previously, to 
a change of identity between r] and p'. 

The meson masses, as function of the density, are plotted in Fig. HUa). The SU(2) 
chiral partners (vr^, a) are now always bound states. The pion is a light quark system for 
all range of densities and the a meson has a strange component, at pp = 0, but becomes 
purely non-strange when Og goes to 35.264°, at pp ~ 3po- At this density the mesons become 
degenerate. This behavior is similar to the non-zero temperature case. 

The SU(2) chiral partner {p,ao) becomes degenerate for 4.0po < Pb < 4.8po, a region 
where they are bound states. In the same range of densities (r/, oq) and (tt^, a) are all 
degenerate. Suddenly the p mass separates from the others becoming a purely strange state. 
This is due to the behavior of 6p that, as already referred, changes sign and goes to 35.264°, 
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at pb — 4.9po- On the other hand, the r]\ that starts as an unbound state and becomes 
bound at pB > S.Opo, turns into a purely light quark system and degenerates with 7r°, a and 
Oo mesons. So, the rj and the p' also change identities. Consequently, contrarily to results 
with temperature, and p' are now degenerate. 

Finally we analyze the behavior of charged mesons with density, plotted in Figs. ^1 (b) 
and [ini(d). The hgure shows that the chiral partners (7r+,aQ*') and (tt ,ag ), panel (b), 
become degenerate for pB — 4po; the chiral partners , k+) and {K~ , k“), panel (d), and 
and , it®), panel (c), do not degenerate in the region of densities considered. We 
notice that, while the results for (tt^ , a^) are affected by the dependence of on density, 
we hnd no substantial differences for the kaonic mesons, whether qb is constant or not. In 
order to understand this, let us remember that for the pion and the Og propagators, the 
dependence on the anomaly enters through the effective coupling {gD)s (see Fig. IHl) so, with 
Qd di. decreasing function of the density, this term will affect less and less the meson masses 
as the density increases. Then, the convergence of the mesons reflects the restoration of the 
Ua( 1) symmetry. Since for kaonic mesons the propagators depend on the anomaly through 
the effective coupling ( 5 'd)„ {{ 90 ) 0 )^ fh® anomaly has little effect on the kaonic masses, as 
the density increases, whether is constant or not, due to the strong decrease of the mass 
of the non-strange quarks. The dominant factor for the calculation of the masses of those 
mesons is the mass of the strange quark, which, although decreasing, remains always very 
high. We can say that the restoration of the axial anomaly does not influence the behavior 
of kaons and of its chiral partners. In addition, we remark that the chiral asymmetry m 
is always different from zero in neutron matter, even for high densities. 

We notice that the convergence between the different chiral partners always occurs at 
densities where the mesons are bound states (see Figs. ITTH and ED, i.e., they are collective 
excitations dehned below the respective qq threshold. 

Case III. In this case we postulate the following dependence for go : gnips) = 
(7D(0)exp[—(ps/pg)^]. The topological susceptibility with this dependence is plotted in Fig. 
El dotted line. From Fig. m a) we see that the density dependence for gnipB) used does 
strengthen the phase transitions, like in the hnite temperature and complete symmetric 
matter cases. The masses of the chiral partners (vr®, og) and {p, a) degenerate at very earlier 
values of the density (ps ~ 2.5po), compared with Case II (where pB — 4pg). Now the 
interval where these four mesons are degenerate is bigger: 2.5pg < Pb ^ 4.8pg. Then the 
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FIG. 10: Density dependence of meson masses and of limits of the Dirac sea continua (dotted lines) 
defining qq thresholds for the mesons. The low-lying solutions are also included. The anomaly 
coupling is a Fermi function (Case II). 

rj mass separates from the others becoming a purely strange state and the (tt^, oq, t]', a) 
mesons become again degenerate in mass. 

In this scenario, y (in Fig. IHl dotted line) goes to zero for ps — 3po- The behavior of the 
mixing angles (Fig. IHl dotted lines) is also qualitatively similar to Case II; 6 s goes to the 
ideal mixing angles for ps > 2.5po and 9p also changes sign, however this happens for lower 
densities, ps — 1-Opo- In panel (b) of Fig. [TUlwe verihed that the degeneracy of and 05 *" 
occurs for pp > 2.5po {pB ^ 4.0po for Case II). In panels (c) and (d) we note a strongest 
decrease of the and k^{bP) masses than in Case II. This is the more relevant effect. 

C. Explicit chiral symmetry breaking without Uyi(l) anomaly 

The absence of mixing effects [gp = 0) in the gap equation for the specihc environment 
now considered induces effects that, although in general qualitatively similar to the previous 
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FIG. 11: Density dependence of meson masses and of limits of the Dirac sea continua (dotted lines) 
defining qq thresholds for the mesons. The low-lying solutions are also included. The anomaly 
coupling is a decreasing exponential (Case III). 

cases studied (finite temperature and symmetric quark matter), have relevant differences: 
(i) it is observed a more significant decrease of the constituent quark mass as compared 
to Mu (the chiral asymmetric parameter plotted in Fig. [7| reflects this behavior); (ii) the 
mass of the strange quark remains constant in the range of densities considered, since there 
are no strange quarks in the medium, due to the fact that Mg > fig (see Eq. 1^ . This 
two facts will have relevant consequences for the mesonic behavior to be discussed in the 
sequel. As it can be seen in Fig. da), and similarly to the previous situations without 
anomaly, tt^ and rj are degenerate in mass, as well as Oq and a and, as the density increases, 
the four mesons become degenerate {pb — 4po)- Some meaningfully differences relatives to 
the other cases with go = ^ appear, however, above pb — 5.5po- Due to the absence of 
the anomaly, there are no mixing effects and the mixing angles have, therefore, always ideal 
values. However, we observe a change of sign of the pseudoscalar angle, 6*p, at that density 
{6p = —54.736°, for pp < 5.5po,9p = 35.264°, for pp > 5.5po) a behavior that seems specihc 
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Pb^Po Pb^P, 


FIG. 12: Density dependence of meson masses and of limits of the Dirac sea continua (dotted 
lines) defining qq thresholds for the mesons. The low-lying solutions are also included. The axial 
anomaly is absent {qd = 0). 


of the type of matter under study. This implies, as usual, that the rj meson, non strange 
up to this density, becomes purely strange afterwards, the opposite happening to r/', that 
changes the role with r] from now, being degenerate with 7r°, oq , cr. A consequence of the 
strange quark mass remaining constant is that the mesons with only a strangeness content 
keep their masses constant As in the Cases I, II, and III, the /o-meson shows no 

tendency to become degenerate with any other meson. 

In panel (b) the Uq are always bound states and we verify that the degeneracy of (oq”, vr“) 
and (ao'",vr+) occurs at different baryonic densities, respectively, pb — 3.5po and ps — 4.1po- 
This may indicate the existence of two separate hrst-order phase transitions in the non- 


strange sector, in agreement with the conclusions of |47|. 

In panels (c) and (d) we notice a strongest decrease of masses as compared to 

those of K^. This is due to a more pronounced decrease of with increasing density. In 
addition, the splitting between charge multiplets of pions and kaons is always manifest as 
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expected. 


VII. SUMMARY AND CONCLUSIONS 

In this work we investigated different patterns of restoration of axial symmetry, in con¬ 
nection with the restoration of chiral symmetry, in a model with explicit breaking of the 
Uyi(l) anomaly. The restoration of axial symmetry at non-zero temperatnre (density) has 
been discussed using two different decreasing functions of temperature (density) for the cou¬ 
pling anomaly, gu: one of them is inspired in lattice results (Case II) for the topological 
susceptibility and the other is a simple exponential function (Case III). This results were 
compared with the case were Qd = constant for all temperatures and densities. 

We verihed that in the last case there is always an amount of Ua( 1) symmetry breaking in 
the particle spectrum even when chiral symmetry restoration in the non-strange sector occurs 
at high temperature (density). To complement the information provided by the effective 
restoration of axial symmetry, the extreme case qd = 0 has also been considered. For a 
more complete understanding of the density effects we considered two different scenarios of 
quark matter: (i) symmetric quark matter; and (ii) neutron matter in /3-equilibrium. So, 
the different patterns of axial symmetry in the vacuum state, with Qd = constant (Case I), 
Case II and Case III, and = 0, have been applied in a hot medium, in symmetric quark 
matter and in neutron matter. 

Since in all cases chiral symmetry is explicitly broken by the presence of non-zero current 
quark mass terms, the chiral symmetry is realized through parity doubling rather than by 
massless quarks. So, the identihcation of chiral partners and the study of its convergence is 
the criterion to study the effective restoration of chiral and axial symmetries. An important 
information is also provided by the mixing angles and we verify that, in the scenario of 
effective restoration of axial symmetry, the mixing angles converge to the situation of ideal 
flavor mixing: (i) the a and f] mesons are pure non-strange qq states, while /o and rj' are 
pure strange ss excitations for symmetric matter and non-zero temperature cases; (ii) the rj 
and rj' change identities for neutron matter case. 

In the conditions of explicit breaking of chiral symmetry (real world) we worked, SU(3) 
symmetry is not exact and, even in the limiting case g^ = 0, the strange sector does 
contribute with signihcant effects even at high temperature (density) as it is visible in the 


34 



behavior of /o and rj [r]') mesons. 

We can conclude that in Cases II (or III) the Ua( 1) symmetry is effectively restored 
above the critical transition temperature of the SU(2) chiral phase transition. But, in 
the region of temperatures (densities) studied we do not observe signs indicating a full 
restoration of U(3)®U(3) symmetry as, for instance, the degeneracy of both qq and /o 
mesons with the pion. In fact, as we work in a real world scenario (m^ = rud << rus), 
we only observe the return to symmetries of the classical QCD Lagrangian in the non- 
strange sector. The dynamics of the system at low temperatures or densities is dominated 
by quantum effects of both chiral and Uyi(l) breaking symmetries. This is manifest in the 
low-lying mesonic spectrum. As the temperature or density increase our model simulates, at 
least phenomenologically, features of the large hadron mass spectrum. In such systems both 
chiral and Ua( 1) symmetries must be restored, which is signaled through a systematical 
appearance of degenerate chiral and axial partners. 

We started with explicit symmetry breaking in the presence of the Ua( 1) anomaly in 
the vacuum state, with the axial symmetry being effectively restored by thermal (density) 
effects. The results are based on a schematic model, however, it includes some of the 
main ingredients for a reliable qualitative description of the high temperature or density 
regime of matter. The anomalous effective interaction vanish under extreme conditions of 
temperature/density as required by asymptotic freedom of QCD. A more realistic approach, 
which includes the enlargement of this behavior to the scalar-pseudoscalar interaction, can 
be done in the framework of a model with finite range form factors. 
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APPENDIX A 

In this Appendix we present some technical details of the model formalism in the vacuum 
state and at finite temperature and chemical potential. 
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1. Propagators and polarization operators for pseudoscalar mesons 


The effective c 
bilinear {qX°'q) 


uar 


c Lagrangian Q have been obtained making a contraction of one 


35, 


36 


37| with the projectors Sab, Pat given by 


Sab QS^ab T doPabc {qX q) , (-^1) 

Pab QS^ab doPabc {qX q) , (-^ 2 ) 


where {qX'^q) are vacuum expectation values. The constants Dabc coincide with the SU(3) 
structure constants dabc for a,b,c = (1, 2,, 8) and D^ab = -Dqoo = 

The effective model Lagrangian Q has been written in a form suitable for the usual 
bosonization procedure. This can be done by the integration over the quark fields in the 
functional integral. So, the natural degrees of freedom of low-energy QCD in the mesonic 
sector are achieved. It gives the following effective action: 




— iTr In 


-m + craA“ (i75)(</5aA“) 


(A3) 


The notation Tr stands for the trace operation over discrete indices {Nj and and 
integration over momentum. The helds cr“ and 93 “ are scalar and pseudoscalar meson nonets, 
respectively. 

To calculate the meson mass spectrum, we expand the effective action (IA3fl over meson 
fields. Keeping the pseudoscalar mesons only, we have the effective meson action 

w-'./lM = - nS.(-P)] = -\‘p'‘(D^AP)rV, (A4) 

with being the polarization operator, which in the momentum space has the form 

nS,(P) = tN, j [Si(p)(A“)«(n5)S,(p + P)(A‘),.(n5)] , (as) 

where iin is the trace over Dirac matrices. The expression in square brackets in (IA4I1 is the 
inverse non-normalized meson propagator (P^(P))“^. 

For the non-diagonal mesons tt , iF, the polarization operator takes the form 


npPo) = 4 ((/; + /;) - - (m, - M,f] 7(p„)), (as) 
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where the integrals I{ and /^■^(Pq) are given by 




= iN, 


d^p 1 

(27r)"^ _ ]\^2 


Nc 

4vr2io P, ’ 


(A7) 


jij(p N .JV f _i_ 

V {2t,Y {p^ - MMp + P,Y - M]) 

_ Nc p^dp Ei + Ej 

Jo P2 - (P, + P,)2 ’ 


(A 8 ) 


where Eij = y p^ + MA is the quark energy. To regularize the integrals we introduce the 
3-diniensional cut-off parameter A. When Pq > Mi + Mj it is necessary to take into account 
the imaginary part of the second integral. It may be found, with help of the ie -prescription 
Pq — Pq — ie. Using 

lim -= V — \-i7rS{y) (A9) 

e^o+ y — ze y 

we obtain the integral 


IHPo) 


^ p'^dp Ei + Ej . p* 

^ P2-(P, + P,)2 {E* + E*y 


(AlO) 


with the momentum: p* = {Pq — {Mi — MjY){Po — {Mi -|- MjY)/2Po and the energy: 

Eh = y'(p-)" + Mh. 

To consider the diagonal mesons tt'^, p and p' we take into account the matrix structure 
of the propagator in (IA4fl . In the basis of tt^ — p — p' system we write the projector Pab and 
the polarization operator as matrices: 



P 33 P 30 P 38 


/ TT^ TT^ TT^ \ 

Pab — 

Po3 Poo P 08 

and = 

TT^ TT^ TT^ 

'*^■'-03 ■‘^'‘-00 '‘-'-08 


y P83 P80 P88 j 


Trn Tin TT^ 
y J^J^83 J^J^80 J^J^88 y 


(All) 


The non-diagonal matrix elements P 30 = -^gD{{ququ) - {qdqd)), P 38 = -^goiiququ) - 
(qdqd)), Uso = V^[nL(Ab) - ^dd{Po)] and Ugg = l/V3[U^^{Po) - n^rf(Po)] correspond to 
71^ — p and 7 r° — p' mixing. In the case {qu qu) = {dd qd), the 7 r° is decoupled from the p — p' 
system and the preceding matrices have the non-vanishing elements: 
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and 


where 


P33 = gs + gn {qs Qs) , 


^00 = gs- ^gn {{qn Qu) + {qd qd) + {qs qs )), 

P33 = gs + ^gn (2 {qu qu) + 2 {qa qd) - {qs qs)), 


P 08 — -Psn — 


80 


3^2 


go {{qu qu) + {qd qd) - 2 {qs qs)) , 


nfo(fo) = I [^.(Fo) + nJ,(F„) + nf,(fi,)], 
n£,(Fo) = I [nL(Po) + nZiPo) + 4nf.(p„)], 

<(p«) = nJ,(F„) = ^ [n^„(p„) + n^(p„) - 2 nf.(p„)], 


n?(p„) = 4(2/; - p|/“(p„)). 


(A12) 

(A13) 

(AM) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 


The procedure to describe scalar mesons is analogous. We present below the most relevant 
steps. 

To calculate the meson mass spectrum, we expand the effective action (El over meson 
fields. Keeping now the scalar mesons only, we have the effective meson action 

ibSM = K - nf,(P)] <t‘ = -l^“(Df,(P))-V‘, (A20) 

with nf^(P) being the polarization operator, which in the momentum space has the form of 
dMI) with (* 75 ) substituted by the identity matrix. 

The polarization operator associated with the non-diagonal mesons (oq, u, /o) has the 
form 

ng(P„) = 4 ((/• + I() + [PS - (M7 + Mf)] /«(P„)) . (A21) 


To consider the diagonal mesons Oq , a and /o we take into account the matrix structure 
of the propagator in (jA20jl . In the basis of — a — fo system we write the projector Sab 
and the polarization operator flf^ as matrices: 


/ q q q \ 
*J33 *->30 *->38 


Sab = 


*5*03 *500 *508 
*583 *580 *588 


nfa nf„ n| ^ 


and nfb = 


/ 


^03 Hqq flgg 

nfs nfo nfg 


(A 22 ) 


/ 
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ab 


and 


ab 


nss 0 

0 nf;, 


In the case g„) = {qd qd) the preceding form of the matrices is rednced to 

533 0 ' ' 

^ 0 Sab 

with 

533 = 9s- 9d {Qs Qs) , 

2 _ 

5oo = 9s + ^9 d {{qu qu) + {qdqd) + qs )), 

588 = 9s- ^9 d (2 {qu qu) + 2 {qd qd) - fe ?*)), 


5o8 — 580 — 


3^/2 


9d {{qu qu) + {qd qd) - 2 (g* g^)) 


Analogously, we get 


nS)(a) = 3 [nfJPo) + nf,(p„) + nf,(F„)], 
n|(a) = I [njja) + nf,(p„) + 4 nf.(p„)], 
nS,(a)) = nf„(F„) = ^ [ni(p„) + n®,(p„) - 2nf,(P(, 


where 


We also obtain 


and 


n,^(Po) = 4(2/i + [p2 _ 4Mf]J“(Po)). 
D-^ = {A + C)- V(C - Ay + 4P2 
P 7 ; = {A + C) + ^/{C-Ay + 4B\ 


(A23) 


(A24) 

(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 

(A32) 

(A33) 


where the expressions for A, B and C are formally analogous to those for pseudoscalars. 

The masses of the a and /o meson can now be determined by the conditions P7^(Mo-, 0) = 
0 and Dj^{Mf^, 0 ) = 0 . 

2. Model formalism at finite temperature and chemical potential 


The NJL model can be generalized to the hnite temperature and chemical potential case. 
It can be done by the substitution Q 

f d'^p 1 f dSp 


(2vr) 


—ip J (27r)3 




(A34) 
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where (3 = 1/T, T is the temperature and the sum is done over Matsubara frequencies 
ujn = (2n + l)7rT, u = 0, ±1, ±2,..so that po —*■ + P with a chemical potential fi. 

Instead of integration over po we have now the sum over Matsubara frequencies which can 
be evaluated 


^ ^ h{uJn) — ^ 

" Rezm^o 


(1 y ) R6s[/i(cu^), 


“b y ( '^m ) RgS [h 


where f{z) and f{z) are the Fermi distribution functions for quarks and antiquarks: 


/W = 


1 


/» = 


(A35) 


(A36) 


X -|_ Q^iz-u) ’ ^ V / g/3(z+tj) ■ 

As 1 — f{z) = f{—z), we introduce, for convenience, the Fermi distribution functions for the 
positive (negative) energy state of the ith quark: 

1 


"i =/i(±Bi) = Y-+j±s(E.w)' 
At hnite temperature the integral II dSl) takes the form 


(A37) 


IliT,fi^) = 


Nr 


P^dp 


- rii ) 


(A38) 


47r2 J Ei 

The integral I 2 {P) depends now on the temperature T and two chemical potentials pi, pj 
which are appropriated to quark flavors 


Il\Po,T, Hi, ij,j) = -Nc 


Pp 

1 


2E,{Ei + Po-{Hi-I^j)y-E] 

1 


n; 


2Ei (Ei — Pq + (hi ~ hj))^ ~ E. 


2 


nt 


2P,(P,-Po + (p,-p,))2-Pf ^ 

1 1 


n. 


2 E,{E,+Po-{Hi-HjW-Ef ^ 
For the case i = j, with imaginary part, we have the expression 


lf{Po,T,Hi) = 


27r2 


V 


p^dp 


Ei Pd-4Ef 


2 V * 


n- — n. 


-) 




(A39) 


(A40) 
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Having these integrals as functions of the temperature and chemical potentials, we can 
investigate the meson properties in hot/dense matter. 

3. Topological susceptibility 

The topological susceptibility is given by 


= j 


d^a; e-*'^"(0|TQ(x)Q(0)|0)_ted, 


(A41) 


where Q{x) is the topological charge density. The general expression for y in NJL model 
has been obtained in m 


X = Wd 2n^„(0) {Qu quf {Qs Qsf + n£(0) {Qu quf 



(A42) 
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